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ABSTRACT 

We investigate the optical emission-line flux ratios of narrow-line regions, in order to determine whether 
the formation of AGN jets requires specific accretion conditions. We find that bright compact radio 
galaxies, which are powerful radio galaxies in the early stage of the jet activity, exhibit systematically 
larger flux ratios of [Oi]A6300/[Oni]A5007 and smaller flux ratios of [Oiii]A5007/[Oin]A4363 than radio- 
quiet (RQ) Seyfert 2 galaxies. Comparing the observed line ratios with photoionization models, it is 
found that the difference in the flux ratio of low- to high- ionization lines (e.g., [Oi]A6300/[Oiii]A5007) 
can be well understood by the difference in the spectral energy distribution (SED) of ionizing sources. 
Powerful young radio-loud (YRL) AGNs favor SED without a strong big blue bump, i.e., a radiatively 
inefficient accretion flow (RIAF), while RQ AGNs are consistent with the models adopting SED with 
a strong big blue bump, i.e., a geometrically thin, optically thick disk. These findings imply that the 
formation of powerful AGN jets requires the accretion disk with harder ionizing SED (i.e., a RIAF). 
We discuss the obscuring structure of YRL AGNs as a plausible origin of the difference in flux ratios of 
[Oiii]A5007/[Oiii]A4363. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies:jets 



1. INTRODUCTION 

The formation of relativistic jets is one of the funda- 
mental issues in active galactic nuclei (AGNs) physics. 
The question how the jets form is relevant to supermassive 
black hole (BH) cases as well as X-ray binaries (XRBs; see 
reviews by Fender 2003). Thanks to their short dynamical 
timcscale, XRBs in various states have been observed in 
detail, and it is found that individual sources occupy par- 
ticular accretion state with various X-ray luminosity (e.g., 
Fender et al. 2004; Remilland & McClintock 2006). These 
observations suggested that at the low luminosity (i.e., the 
low mass accretion rate) the jet power is a monotonic func- 
tion of the accretion luminosity. On the other hand, the 
jet power may change abruptly by a large factor, following 
transitions between the soft state (a spectra characterized 
by the black body components) and hard state (a spectra 
characterized by the power-law components) at the high 
luminosity (i.e., the high mass accretion rate). By analogy 
to XRBs, it is expected that the power of AGN jets also 
depends on the mass accretion rates in a similar way. For 
example, Rawlings & Saunders (1991) found that for the 
extended radio galaxies whose projected linear size (LS) is 
larger than 10 kpc the jet power is proportional to the opti- 
cal narrow line luminosity of NLRs which are photoionized 
by the nuclear continuum radiation (see also Willott et al. 
1999). Ghisellini & Celotti (2001) claimed that the less 
powerful radio galaxies such as Fanaroff- Riley type I radio 
galaxies (FR Is) exhibit the low accretion rate, while the 
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powerful radio galaxies such as Fanaroff-Riley type II radio 
galaxies (FR lis) show the high accretion rate. This work is 
consistent with a scenario (e.g., Marchesini, Celotti, & Fer- 
rarese 2004; Wold et al. 2007; Hardcastle et al. 2007) that 
FR Is have a radiatively inefficient accretion disk (RIAF) 
(e.g., Narayan & Yi 1995), while FR Us have a standard 
accretion disk (e.g., Shakura & Sunyaev 1973). 

In investigating the disk-jet connection in AGNs, how- 
ever, we should keep in mind the relatively longer life time 
of the extended radio sources, i.e., FR Is and FR Us. The 
reasons are as follows. The time-scale of the transition 
to the different physical states of the accretion disk is 
At(XRBs) w 10~ 2 - 10 _1 yr (a few days - a few weeks 
change) for XRBs (e.g., Fender et al. 2004). Assuming 
Mbh(XRBs) = 10M Q and Mbh (AGNs) = 1O 8 M , the 
transition time-scale for AGNs is Ai(AGNs) « 10 5 ~ 6 yr 
because of At oc Mbh- Thus, the physical states of the 
accretion disk may change in less than the typical age of 
extended radio galaxies, i.e., t agc rj 10 7 ~ 8 yr (e.g., Parma 
et al. 1999; Bird, Martini, & Kaiser 2008). In other words, 
the observed relation between the jet power and the ac- 
cretion power for the extended radio galaxies may not be 
proving the direct disk-jet connection in AGNs. Thus, in 
order to investigate the disk-jet connection in AGNs more 
directly, it is essential to examine the physical states of 
the accretion disk in younger radio galaxies. 

The young and compact radio galaxies such as compact 
symmetric objects (CSOs; LS < lkpc) and medium-size 
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symmetric objects (MSOs; LS =1 — 10 kpc) have been 
discovered by several authors (e.g., Wikinson ct al. 1994; 
Fanti et al. 1995; Readhead et al. 1996). Among these 
compact radio sources, the bright objects show the FR I 
I- like morphology along with hot spots (reverse shocks), 
which are considered to be a signature of supersonic ex- 
pansions (Owsianik, Conway & Polatidis 1998; Taylor et 
al. 2000; Tschager et al. 2000; Giroletti et al. 2003; Po- 
latidis & Conway 2003; Gugliucci et al. 2005; Nagai et al. 
2006; Gugliucci et al. 2007; Luo et al. 2007). In contrast, 
the faint compact radio sources do not show hot spots (e.g., 
Kunert-Bajraszewskaet al. 2005; Giroletti 2007) and some 
of them show the FR I- like morphology (Giroletti 2007). 
Since the age of CSOs and MSOs is less than ~ 10 5 yr, they 
are recognized as newly born young AGN jets (e.g., Fanti 
et al. 1995; Readhead et al. 1996, O'Dea & Baum 1997; 
Stanghcllini et al. 1998; Snellen ct al. 2000; Dallacasa et 
al. 2000; Orienti et al. 2007). Thus, CSOs and MSOs 
are adequate targets to examine the disk-jet connection 
in AGNs since the time lag between the optical and ra- 
dio activity is relatively small (Ai(AGNs) < 10 5-6 yr). In 
this paper, we will use bright CSOs and MSOs (hereafter 
powerful young radio-loud [YRL] AGNs) to investigate the 
disk-jet connection in AGNs. 

In order to probe whether the formation of AGN jets 
requires any specific accretion conditions, it is essential to 
determine the physical states of the accretion disk of pow- 
erful YRL AGNs. To this end, it can be useful to examine 
a systematic difference in the SEDs produced by accretion 
disks of powerful YRL AGNs and radio-quiet (RQ) AGNs. 
However, it is challenging to derive the intrinsic SEDs of 
powerful YRL AGNs from the hard X-ray observations, 
because of the contamination of AGN jets and the obscura- 
tion of the accretion disk by the dusty torus. Alternatively, 
the optical emission-line diagnostics of narrow-line region 
(NLR) can be used to distinguish the physical states of the 
accretion disk, since the NLR is thought to be photoion- 
ized by the nuclear radiation for RQ Seyfert galaxies (e.g., 
Yee 1980; Shuder & Osterbrock 1981; Evans et al. 1999) 
and even for extended radio galaxies (e.g., Villar-Martin 
et al. 1997; Nagao et al. 2006). However, very little is 
currently known about the optical emission-line proper- 
ties of powerful YRL AGNs. By examining the difference 
of optical line ratios between powerful YRL AGNs and the 
extended radio galaxies, Morganti et al. (1997) found that 
powerful YRL AGNs show the characteristics very similar 
to that of extended radio galaxies. However, it is not clear 
how the physical state of accretion disks (especially SEDs 
of accretion disks) is different between the powerful YRL 
AGNs and the RQ AGNs. Answering this question will 
provide crucial information on the disk-jet connection in 
AGNs. 

The main goal of this paper is to investigate whether 
there is the difference in SEDs between the powerful YRL 
AGNs and Seyfert 2 galaxies (i.e., RQ AGNs), by utiliz- 
ing the optical emission-line diagnostics of NLRs. We col- 
lected samples of powerful YRL AGNs and RQ AGNs (§2), 
and generated various photoionization models with differ- 
ent input SEDs (§3). In §4, we compared the observational 
data with models. On the basis of our findings, Discussion 
on accretion conditions for the formation of powerful AGN 
jets is presented in §5, and Conclusions follow in §6. 



2. DATA COMPILATION 

We focus on the properties of ionized gas in NLRs for 
powerful YRL AGNs to investigate the disk-jet connec- 
tion. To reduce possible inclination effects, it is important 
to compare AGNs with the same optical type, since the 
optical narrow emission line ratios are often affected by 
the geometrical effects (e.g., Murayama & Taniguchi 1998; 
Schmitt 1998; Nagao, Murayama, & Taniguchi 2000; Na- 
gao, Murayama, & Taniguchi 2001a). Out of 219 radio 
galaxies available in the literature, we compiled 32 YRL 
AGNs (20 narrow-line AGNs and 12 broad-line AGNs) 
with the optical spectra, from which the stellar contin- 
uum is subtracted, and high radio luminosity, logP5GHz > 
25.0 W/Hz, using papers by Pearson & Readhead (1988), 
Tadhunter et al. (1993) and O'Dea (1998). Since only few 
broad-line YRL AGNs have the detected emission lines 
(e.g., [Oi]A6300 line, [On]A3727 and [On]A5007) which are 
essential for our study, we finally selected 20 narrow-line 
YRL AGNs with these lines. In other words, the YRL 
AGN sample consists only of radio galaxies. The prop- 
erties of powerful YRL AGNs arc presented in Table 1. 
The mean projected linear size of selected radio sources 
is 3.2 kpc, which is small enough to discuss the disk-jets 
connection in AGNs (see §1). The average radio power at 
5 GHz, log P5GHZ is ~ 27.0 W/Hz, comparable to power- 
ful extended radio galaxies (Laing, Riley & Longair 1983). 
Note that the radio power of the faint YRL AGNs, ex- 
cluded from our sample, is logPsGHz < 25.0 W/Hz (e.g., 
Giroletti 2007). We use only the forbidden lines to avoid 
the difficulty of separating the broad component and the 
narrow component of the recombination lines (e.g., Ha and 
H/3 lines). Then, we compiled the flux ratios of forbidden 
emission- lines from the literature (Table 2). The typical 
errors of each line are « 10 — 30%. 

To compare with narrow-line YRL AGNs, we selected 
a sample of RQ Seyfert 2 galaxies from the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) archive, using the 
emission- line flux catalog made by the MPA/JHU group 
(Kauffmann et al. 2003). We selected 624 RQ Seyferts, 
of which First Images of the Radio Sky at Twenty cm 
(FIRST) flux is lower than 1 mJy and optical emission- 
line fluxes are well- measured from high S/N (> 10) SDSS 
spectra. 

In this paper, we do not evaluate the amount of ex- 
tinction of host galaxies based on the Balmer decrement 
method (e.g., Osterbrock 1989), since for only few sources 
in the powerful YRL AGN sample detected were both the 
narrow Ha and H/3 lines. The data presented in Table 
2 are not corrected for Galactic extinction. The effect of 
dust extinction on our results will be discussed in §4.1. 

3. PHOTOIONIZATION MODEL 

To characterize SEDs produced by the accretion disk 
in powerful YRL AGNs, we will compare emission line 
properties with photoionization models. We will examine 
whether there is a systematic difference in SEDs from ac- 
cretion disks between the powerful YRL AGNs and the RQ 
AGNs. The continuum radiation from the accretion disk 
photoionizes the NLR, hence, the difference in the SED 
can be related to the difference in the physical states of the 
NLR. Then, we compare the observed flux ratios with pho- 
toionization model calculations, using the publicly avail- 
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able code CLOUDY, version 06.02c (Ferland 1998). 

The method of our calculations is based on solving the 
equations of statistical equilibrium, the ionization and 
thermal equilibrium through gas clouds as a function of 
depth. Then, we obtain the ionization structure of gas in 
NLRs and the volume emissivity of the forbidden emis- 
sion lines (i.e., [Oi]A6300, [On]A3727, [Sn]A6717+6731, 
[Om]A 5007 and [Neili]A3869). Here we briefly describe 
the adopted parameterizations. We assume uniform den- 
sity gas clouds with a plane-parallel geometry. The pa- 
rameters for the calculations are (i) the hydrogen density 
(nn); (ii) the ionization parameter (U = Q(H)/47rcr 2 nH, 
where Q(H) is the number of the ionizing photons emitted 
by AGNs per second), which is defined as the ratio of the 
ionizing density to the hydrogen density; (iii) the chemical 
composition of the gas; (iv) the shape of SEDs of the input 
continuum radiation. While the gas metallicity is assumed 
to be solar (Grevesse & Sauval 1998; Holweger 2001), we 
carry out a series of model runs covering various density 
and ionization parameters: 10 2 cm~ 3 < njj < 10 5 cm~ 3 
and 10- 4 < U < lO" 15 . 

For the input continuum radiation, we consider two 
types of SEDs. One is a typical SED produced by a stan- 
dard accretion disk (Shakura & Sunyaev 1973), which is 
characterized by a strong big blue bump (BBB) in the 
wavelength range of UV to soft X-ray. As a template of 
this type of SED (hereafter, SED with BBB), we adopt 
the empirical expression using the following function (see 
Ferland 1997; Nagao et al. 2001b): 

where the infrared cutoff of the BBB component, fcTiR = 
0.01 ryd (lryd = 13.6 eV), the characteristic temperature 
of the BBB, T B b = 4.9 x 10 5 K, the slope of the low-energy 
side of the BBB, «uv = —0.5, the UV to X-ray spectral 
index, a ox = —1.35 and the slope of the X-ray power- 
law continuum, a x = —0.85 which are typical values of 
Seyfert galaxies (see Nagao et al. 2001b for more details). 
Since the characteristic temperature of the BBB in some 
extreme AGNs, e.g., narrow- line Seyfert 1 galaxies, can 
reach at Tbb ~ 10 6 K, we will discuss the effects of chang- 
ing Tbb on the line flux ratios. Note that the parameter 
a in cq. (1) is determined from the adopted value of a ox . 
The last term in equation (1) is set to zero below 1.36 eV, 
and the continuum is assumed to fall off as above 100 
keV. 

The other kind of SED is a typical SED produced by 
a RIAF, which is expected to show a power-law with no 
strong BBB as shown by various works (e.g., Narayan & Yi 
1995; Chen, Abramowicz, & Lasota 1997; Narayan, Kato, 
& Honma 1997; Manmoto, Mineshige & Kusunose 1997; 
Kurpiewski & Jaroszynski 1999; Quataert & Narayan 
1999; Manmoto 2000; Kino, Kaburaki, & Yamazaki 2000; 
Ohsuga, Kato & Mineshige 2005; Yuan, Ma & Narayan 
2008). We adopt here the SED model presented by 
Kurpiewski & Jaroszynski (1999) as the template for SEDs 
of RIAF (hereafter SED without BBB). This SED can be 
expressed by a single power-law continuum in the range 
of 10 12 Hz to 10 20 Hz, with the exponential cutoffs at 
10~ 4 ryd and 10 4 ryd. We choose the photon index of 
«pl = 0.89 (f v = ^~ apL ), as predicted for the case of a 
non-rotating black hole. Since there are various shapes 



of the calculated SED for a RIAF (i.e., various opl), we 
will check the effect of altering q;pl in the wide range of 
0.78 < a PL < 1.0. 

Using these two SED templates for a standard disk and 
RIAF, respectively (see Figure 3 in Nagao et al. 2002), we 
calculate emission line flux ratios using CLOUDY. These 
simple SED templates can be used to elucidate the effect 
of the presence or absence of the BBB component, which 
is the main difference in the SED between a standard disk 
and RIAF. Note that these template might be too simple 
to generate accurate narrow-emission line fluxes in AGNs. 
The calculations are stopped when the temperature falls 
below 3000 K, where the gas docs not contribute signifi- 
cantly to the observed optical emission line spectra. 

4. RESULTS 

In this section, we compare the photoionization model 
calculations with observed emission line properties of pow- 
erful YRL AGNs and Seyfert 2 galaxies. 

4.1. Comparison with models with/without BBB 

The emission line ratio, [Oi]A6300/[Oiii]A5007, of the 
NLR can be used to characterize the SED of accretion 
disk since the line ratio depends on the shape of SED 
of the photoionizing source. The harder spectra (with- 
out BBB) of optically thin disks create larger ratios 
[Ol]A6300/[Om]A5007 due to the following reasons. The 
mean free path of ionizing photons at higher energies be- 
comes longer, since the cross section for ionization of hy- 
drogen (er(H)) is given by a (H) cx (u/v^)^ for v > j^l, 
where v\, is the frequency of the Lyman limit. With a flat 
power-law spectrum (apL < 1), a relatively large num- 
ber of photons are available at higher energies even after 
the photons near the Lyman limit are all absorbed around 
the ionization front. Then, the higher energy photons can 
reach at larger radii from the ionization source, hence, the 
extended partially ionized region can be formed. Further- 
more, the low-ionization emission arises in the extended 
partially ionized region rather than the fully ionized re- 
gion where the abundance of neutral oxygen is negligible. 
Therefore, the volume emissivity of low-ionization emis- 
sion lines (i.e., [Oi] A6300) can be higher at larger radii from 
the ionization sources, compared to that of high-ionization 
emission lines (i.e., [Om]A5007) (see Figure 5 of Nagao et 
al. 2002). Note that the ionization potential of neutral 
oxygen is nearly identical to that of hydrogen. 

Figure 1 represents emission line ratios ([Oi]A6300/[Oni] 
A5007 versus [On] A3727/[Olll]A5007) from our model 
calculations and observations. The horizontal axis 
([Oii]A3727/[Oiii]A5007) indicates the difference of n H be- 
cause the critical density of [On]A3727 (rin = 4.5 x 
10 3 cm~ 3 ) is far smaller than that of [Oiii]A5007 (nn = 
7.0 x 10 5 cm -3 ). This set is free from the chemical abun- 
dance effect. As shown in Figure 1, powerful YRL AGNs 
have relatively higher [Oi]A6300/[Oiii]A5007 than Seyfert 
2 galaxies, implying that there is a difference in the shape 
of ionizing continuum. The observed line ratios were not 
corrected for both the Galactic extinction or internal ex- 
tinction of host galaxies. Instead, we indicate by the arrow 
(upper left in Figure 1) the effect of dust extinction when 
the dust-extinction correction of A v = 1.0 mag is applied, 
assuming a standard Galactic reddening law (Cardelli et 
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al. 1989). For type 1 AGNs, the dust reddening of NLR 
emission is small (A v < 1.0 mag; e.g., Rodriguez- Ardila 
et al. 2000). The dust-extinction of NLRs for type 2 
AGNs is approximately 1.0 mag larger than that for type 1 
AGNs (e.g., De Zotti & Gaskell 1985; Dahari & De Rober- 
tis 1988). 

The models using SED without BBB predict higher 
[Oi]A6300/[Oiii]A5007 ratio than the models using SED 
with BBB. The trend of the model prediction coincides 
with the observations. The line ratios of Seyfert 2 galaxies 
are consistent with models adopting the SED with BBB 
in the range of 10 -35 < U < 10~ 30 and 10 40 cm _3 < 
«h < 10 5 '°cm~ 3 , while the line ratios of powerful YRL 
AGNs can be explained by the model adopting the SED 
without BBB. Note that the derived hr may not reflect 
the actual number density of narrow-line clouds because 
of the dust extinction effects. However, the direction of 
extinction correction is perpendicular to the direction of 
the difference between powerful YRL AGNs and Seyfert 
2 galaxies, indicating that the observed difference is not 
due to the difference in the degree of the dust extinction 
between powerful YRL AGNs and Seyfert 2 galaxies. 

In Figure 2, we present the frequency distributions of 
log([Oi]A6300/[Oin]A5007) for the powerful YRL AGNs 
(blue) and Seyfert 2 galaxies (red). The mean value 
of log([Oi]A6300/[Om]A5007) for powerful YRL AGNs is 
-0.36 ±0.26, while that of log([Oi]A6300/[Om]A5007) for 
Seyfert 2 galaxies is —1.17 ± 0.05. The Kolmogorov- 
Smirnov (K-S) test indicates 10 -4 probability that the 
observed frequency distributions of [Oi]A6300/[Oiii]A5007 
of powerful YRL AGNs and Seyfert 2 galaxies originate 
from the same underlying population, suggesting that the 
emission line ratios of [Oi]A6300/[Oiii]A5007 for powerful 
YRL AGNs are statistically larger than those of Seyfert 2 
galaxies. 

In figure 3, we compare the observations 
with the model predictions in other diagnos- 
tic diagrams, i.e., [Oi]A6300/[Oiii]A5007 versus 
[Sii]A6717+6731/[Oiii]A5007 and [Oi]A6300/[Oiii]A5007 
versus [Ne m]A3869/[Oii]A3727. The physical meaning of 
these diagrams is same as Figure 1. Similar to the trend 
in Figure 1, powerful YRL AGNs and Seyfert 2 galaxies 
are well separated in these diagnostic diagrams, suggest- 
ing that the different line ratios of powerful YRL AGNs 
can be explained by models with the photoionizing SED 
without BBB. However, one may doubt that the difference 
of SED is not a unique solution to explain the observed 
differences. Thus, we discuss the other possibilities for 
high [Oi] A6300/[Om]A5007 ratios in §5.1. 

4.2. Dependence on the model parameters 

Here, we examine the dependence of emission line ra- 
tios on the model input parameters, Tbb and apL- Fig- 
ure 4 shows the optical emission line ratios ([Oi]A6300/ 
[Olll]A5007 and [Oii]A3727/[Oiii] A5007) as a function of 
apL (left panel) and Tbb (right panel), assuming rtu = 
10 4 cm~ 3 . We find that at given U the variations of Tbb 
and apL does not change significantly the line ratios of 
[Ol]A6300/[Om]A5007, compared to the effect of the dif- 
ference in SEDs of ionizing radiation. Since the change of 
emission line ratios due to «pl or Tbb is relatively small, 



the effect of changing Tbb and apL cannot account for 
the observed difference between powerful YRL AGNs and 
Seyfert 2 galaxies. The effect of both Tbb and a does not 
also change the results shown in Figure 3. Therefore, the 
difference in the observed emission line ratios in Figure 
1 and Figure 3 is consistent with the idea that the NLR 
of powerful YRL AGNs is photoionized by SEDs without 
BBB, while that of the RQ AGNs is photoionized by SEDs 
with BBB. 

4.3. [Om]X5007/[Oui]X4363 

In Figure 5, we investigate [Om]A5007/[Olll]A4363 line 
ratios, which is one of the gas temperature indicators (Os- 
terbrock 1989). The line ratio ([Oiii]A5007/[Oiii]A4363) 
of powerful YRL AGNs is smaller than that of Seyfert 2 
galaxies. The mean value of log([Oin]A5007/[Oiil]A4363) 
for powerful YRL AGNs is 1.50 ± 0.29, while that for 
Seyfert 2 galaxies is 2.03 ± 0.16. The frequency distri- 
butions of these line ratios for the powerful YRL AGNs 
(blue) and Seyfert 2 galaxies (red) are presented in Figure 
6. The K-S test indicates 10~ 4 probability that the ob- 
served frequency distributions of [Olll]A5007/[Oni]A4363 
of YRL AGNs and Seyfert 2 galaxies originate in the same 
underlying population. We will discuss potential physi- 
cal reasons for smaller [Oiii]A5007/[Om]A4363 of powerful 
YRL AGNs in §5.3. 

5. DISCUSSION 

5.1. Other possibilities for high [Oi]X6300/[Om]X5007 

ratios 

We here consider three other plausible possibilities 
which may explain the observed difference in line ratios. 

i) Shock ionization. We consider the shock ionization 
due to the interaction between AGN jets and the clouds 
in NLRs. The strong low-ionization emission lines such 
[Oi]A6300 causes by the shock- heated gas (e.g., Mouri, 
Kawara, & Taniguchi 2000 and references therein) be- 
cause shocks generate a large partially ionized region in 
the gas. To investigate whether or not the shock heat- 
ing by AGN jets is responsible for the difference in the 
narrow emission-line ratio between powerful YRL AGNs 
and RQ AGNs, we compare the observed line ratios with 
shock models in Figure 7. The adopted shock models 
with/ without considering the effect of precursor are pre- 
sented by Dopita & Sutherland (1995). The shock models 
predict too small flux ratio of [Oi]A6300/[Oiii]A5007, too 
small flux ratio of [Nelll]A3869/[On]A3727 and too large 
[Oii]A3727/[Oiii]A5007, compared with the observed line 
ratios of powerful YRL AGNs. To explain all observations 
(Figure 1 and Figure 3) by shock models, a fine-tunning of 
model parameters should be required. Thus, these results 
suggest that the difference in the narrow-line flux ratios of 
powerful YRL AGNs are not mainly caused by the shock- 
heated gas in NLRs. However, we should mention that 
for one of YRL AGNs in our sample (3C 303.1) the shock 
ionization cannot be negligible (e.g., Labiano 2005). Note 
that the NLR properties of powerful YRL AGNs can be 
explained by the photoionization model (on average), in 
spite of the signature of the strong kinematic interaction 
with surrounding gas (e.g., Axon et al. 2000; O'Dea et al. 
2002; Holt et al. 2008; Labiano et al. 2008). 



This effect is that the unshocked gas is ionized by the radiation emitted by the hot shocked gas. 
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ii) Gas density. The systematic difference in the gas 
density in NLR may cause the different line ratios, be- 
cause the critical density of [Oi]A6300 transition (tir = 
1.8 x 10 6 cm~ 3 ) is larger than that of [Oiii]A5007 transi- 
tion (nn = 7.0 x 10 5 cm" 3 ). This possibility can be easily 
rejected as shown in Figure 1 and Figure 3. At given ion- 
ization parameter U, the ratio of [Oi]A6300/[Oiii]A5007 
does not depend on the gas density, if we employ only 
the models with SEDs with BBB. This is because both 
critical densities of [Oi]A6300 transition and [Oiii]A5007 
transition are higher than the typical gas density of NLRs 
in AGNs (10 2 cm- 3 < n H < 10 5 cm- 3 ). This indi- 
cates that the observed difference in the flux ratio of 
[Oi]A6300/[Oiii]A5007 is not caused only by a systematic 
difference in the gas density of NLRs. 

iii) Ionization parameter. We consider the difference in 
the ionization parameter U between powerful YRL AGNs 
and RQ AGNs. Since the effect of the SED difference and 
changing the ionization parameter (i.e., the ratio of the 
ionizing photon density to the hydrogen density) are often 
degenerated (see Nagao et al. 2002), a systematic differ- 
ence in the ionization parameter of gas in NLRs between 
powerful YRL AGNs and RQ AGNs can explain the dif- 
ference in the emission line ratios [Oi]A6300/[Oiii]A5007 
when we fix «h- As shown in Figure 1 and Figure 3 (bot- 
tom), we find that only different U cannot explain the 
difference of line ratios between powerful YRL AGNs and 
RQ AGNs. However, it is possible that the combination of 
riH and U can reproduce the difference between them even 
if we assume the SEDs with BBB. If the NLR of power- 
ful YRL AGNs has generally higher tir and lower U than 
that of Seyfert 2 galaxies, the observed difference can be 
explained without invoking different shape of photoioniz- 
ing SEDs. This possibility cannot be completely ruled out 
from Figure 1 and Figure 3, however, a special fine-tunning 
of model parameters is required to generate the observed 
line ratios of powerful YRL AGNs. 

The degeneracy of the effect of SED shape and of 
the ionization parameter is a limitation in interpreting 
these results. For the future studies, the diagnostic dia- 
gram using the flux ratio of [Arlll]A7136/[Olli]A5007 (e.g., 
[Anii]A7136/[Oni]A5007 versus [Oii]A3727/[Om]A5007 
and [Arm] A7136/[Om]A5007 versus [Neiii]A3727/[Oin] A5007) 
can be used to determine the ionization parameter, since 
[Arlll]A7136/[Olll]A5007 is almost independent of the 
shape of SED in the range of U < 10~ 25 (Nagao et 
al. 2002). Thus, it can be better constrained whether the 
SEDs of YRL AGNs are different from those of Seyfert 
2 galaxies. Unfortunately, these lines are not currently 
available for our sample of YRL AGNs. 

5.2. Disk- jet connection 

In §4, we have shown that the powerful YRL AGNs favor 
SEDs without a strong BBB such as the one expected for 
RIAFs. On the other hand, it is well-known that RIAFs 
are radiativcly inefficient and generate relatively lower Ed- 
dington ratios (e.g., Narayna & Yi 1995; Kato, Fukue & 
Mineshige 1998). 

The adovection-dominated solution exists when the 
mass accretion rate is relatively low as M < M max . If 
we assume a fiducial efficiency is sa 10% (Frank, King 

4 log(i2-i0kcv/i[Oi,i]) = 1-5 ± 0.5 and log(L bol /L 2 _i k e v) = 0.9 ± 0.4 



& Raine 1992), the maximum accretion rate is given by 
-^max ~ 3a 2 LEdd/c 2 at r < 100r g , where r g is the 
Schwarzshild radius of the central black hole, a is the a- 
viscosity (Shakura & Sunyaev 1973) and ^Edd is the Ed- 
dington luminosity corresponding to the Thomson scat- 
tering. The energy release from a RIAF (I/riaf) can 
be obtained as Lriaf ~ 0.1(M/M max )Mc 2 (e.g., Kato, 
Fukue & Mineshige 1998). According to three-dimensional 
magneto-hydrodynamic simulations, a is a < 0.1 (e.g., 
Balbus & Hawley 1991; Machida et al. 2000; Machida & 
Matsumoto 2003). We here adopt the conservative value 
as a = 0.1 and then the maximum luminosity of a RIAF 
at M = M max is given as LRiAF,max ~ 10~ 3 L E dd- 

Thus, it is necessary to consider whether RIAFs are 
able to explain the relatively high [Oiii]A5007 luminos- 
ity (hence, high bolometric luminosity) of powerful YRL 
AGNs. First, we estimated the bolometric luminosity 
of powerful YRL AGNs based on the [Oiii] luminosity. 
The [Oiii] luminosity of 16 powerful YRL AGNs were col- 
lected from Gelderman & Whittle 1994 and Lawrence et 
al. _ 1996 (see Table 1). The average [Oiii] luminosity 
is L [0m ] = 10 42 24±1 - 25 erg/s. To evaluate the bolomet- 
ric luminosity, we adopt the bolometric correction of low 
luminosity AGNs (LLAGNs) whose SEDs can be well ex- 
pressed by a RIAF model because of the lack of BBB (e.g., 
Ho 1999), i.e., Lboi ~ 250L [OlII] with variance of 0.9 dex 
(Panesa et al. 2006) Q. Then, the average bolometric lu- 
minosity of YRL AGNs, Lboi.obs ~ 10 45 erg/s, is required 
to produce the observed [Oiii] line luminosity. This sug- 
gests that Lboi,obs is ~ 10~ 2 L^dd for the massive BHs 
(A/bh = 10 9 Af Q ). Hence, it seems difficult to explain the 
bolometric luminosity of powerful YRL AGNs with the 
RIAF model, because of LRiAF.max O.lLboi.obs- The 
combination of emission line ratios and the bolometric lu- 
minosity estimates implies that the physical nature of the 
accretion disk of YRL AGNs cannot be explained by either 
the classical steady disk, i.e., a standard disk or a RIAF. It 
might imply that the relativistic jet occur in a non steady 
accretion disk. 

Interestingly, the physical states of accretion disks in 
some XRBs (e.g., GX 339—4) associated with relativis- 
tic jets are similar to those of YRL AGNs (Homan et al. 
2005; Remillard 2005), suggesting the similarity of disk-jet 
connection in XRBs and AGNs. This similarity may in- 
dicate that the formation of relativistic jets requires lumi- 
nous accretion disk with harder ionizing SED. However, we 
note that the bolometric luminosity estimated from [OIII] 
luminosity is highly uncertain for both luminous AGNs 
and LLAGNs without BBB. Thus, it is crucial to reduce 
the uncertainty of bolometric luminosity by evaluating the 
bolometric correction using a large sample of broad-line 
YRL AGNs, which are left in our future works. 

5.3. Obscuring structure in powerful YRL AGNs 

We discuss here why the line ratios, logQOin] A5007/[Oiii] 
A4363), of powerful YRL AGNs is smaller than that of 
Seyfert 2 galaxies. As we investigated in §4, several op- 
tical narrow line emission ratios ([Oi]A6300/[Oiii]A5007, 
[Oii]A3727/[Oiii]A5007, [Sii]A6717+6731/[Oiii]A5007 and 
[Nem]A3869/[Oii]A3727) can be well explained by single- 
zone photoionization models. However, it has been well- 
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known that any one-zone models overpredicts the flux 
ratio of [Om] A5007/[Oiii]A4363 (e.g., Koski & Oster- 
brock 1976; Heckman 1980; Ferland & Netzer 1983; Rose 
& Cecil 1983; Keel & Miller 1983; Rose & Tripicco 1984). 

By taking high-density gas clouds into account in 
NLRs, Nagao et al. (2001c) reported that photoioniza- 
tion models can explain low emission-line flux ratios of 
[Oiii]A5007/[Oiii]A4363 - 10, without including shock ef- 
fects. This idea is attributed to the fact that the critical 
density of [Om]A4363 (n H = 3.3 x 10 7 c m - 3 ) is higher 
than that of [Om]A5007 (n B = 7.0 x 10 5 c m - 3 ). In ad- 
dition, they suggested that the [Oiii]A4363 originates in 
the dense gas clouds which are located at the inner re- 
gions compared with [Oiii]A5007, hence, a large fraction 
of gas clouds emitting [Oiii]A4363 line are obscured in 
type 2 AGNs. This idea is supported by the observations 
that the line ratios of [Om]A5007/[Oiii]A4363 of Seyfert 2 
galaxies are larger than those of Seyfert 1 galaxies (Na- 
gao et al. 2001c). If this is the case, the difference in 
[Oiii]A5007/[Oiii]A4363 may reflect the difference in the 
structure of duty torus, i.e., the difference of the geomet- 
rical thickness. We should emphasize that the difference 
in Figure 5 (also Figure 6) is not related to the inclina- 
tion effect, because we compare the same optical types of 
AGNs, i.e., powerful narrow- line YRL AGNs and Seyfert 
2 galaxies. Since the ratio [Oin]A5007/[Olli]A4363 of pow- 
erful YRL AGNs is lower than that of Seyfert 2 galaxies, 
the obscuring torus of powerful YRL AGNs may be thinner 
than that of Seyfert 2 galaxies (see Figure 8). According to 
a recent coevolution model of the BH growth and the cir- 
cumnuclear disk (Kawakatu & Wada 2008; see also Wada 
& Norman 2002), the mass accretion rate onto a central 
BH driven by the turbulent viscosity is higher as the scale 
height of circumnuclear disk increases. This is because the 
angular momentum transfer due to the turbulence is pro- 
portional to the the scale height. To test our scenario, it is 
crucial to examine the difference in the structure of obscur- 
ing torus, using Atacama Large Millimeter/submillimeter 
Array (ALMA). 

In summary, our findings indicate that the powerful 
YRL AGNs favor the accretion disk without BBB, sur- 
rounded by the geometrically thin obscuring torus, while 
Seyfert 2 galaxies prefer the accretion disk with BBB, sur- 
rounded by the geometrically thick torus (see Figure 8). 
These results suggest that the formation of the powerful 
AGN jet is related to the physical states of accretion disks 
and structure of obscuring materials around the accretion 



disk. 

6. CONCLUSIONS 

We examine the optical narrow emission-line flux ratios 
of NLRs, in order to determine whether the formation of 
powerful AGN jets requires specific accretion conditions. 
A sample of powerful YRL AGNs, which are at the early 
stage of the jet activity, and radio-quiet Seyfert 2 galaxies 
are selected. We summarize our main results. 

1. Powerful YRL AGNs exhibit larger flux ratios of 
[Oi]A6300/[Om]A5007 than RQ Seyfert 2 galaxies. 
By comparing the observed line ratios with the pho- 
toionization model calculations, we find that the 
difference between powerful YRL AGNs and Seyfert 
2 galaxies can be well understood by the difference 
in SEDs of ionizing radiation from the accretion 
disk. The powerful YRL AGNs favor SED with- 
out a strong BBB such as the ADAF, while the line 
ratios of Seyfert 2 galaxies are consistent with the 
models adopting SED with a strong BBB (a geo- 
metrically thin, optically thick disk). In contrast, 
the observed difference in the line ratios is difficult 
to explain by the difference in the contribution of 
shocks caused by the AGN jets. 

2. The line ratio [Olll]A5007/[Oni]A4363 of powerful 
YRL AGNs is systematically smaller than that of 
RQ AGNs. If we adopt a scenario that [Oiii]A4363 
originates in the dense gas clouds, which are located 
at the inner regions compared with [Oiii]A5007, it 
is inferred that the obscuring torus around power- 
ful YRL AGNs would be thinner than that around 
Seyfert galaxies. 
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Table 1. Properties of powerful young radio- loud AGNs 



Name 


z 


log P 5 GHz (W/Hz) 


Size (kpc) 


logL [0 iii] (erg/s) 


Hctcrcnccs 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


0023-26 


0.322 


27.0 


4.06 




1 


0316+413 (3C 84) 


0.017 


25.1 


4.6 x 10~ 3 


39.9 


3,5 


0345+337 (3C 93.1) 


0.244 


26.0 


4.06 


40.6 


2,3 


0404+769 (4C 76.03) 


0.598 


27.4 


0.86 


41.9 


2,4 


0518+165 (3C 138) 


0.760 


27.8 


3.94 


44.1 


2 


0605+480 (3C 153) 


0.277 


25.7 


13.1 


41.7 


3,5 


0710+439 


0.518 


27.1 


0.14 


42.4 


2,4 


0954+658 


0.899 


26.4 


~ 1 


43.1 


4,6 


1031+567 


0.459 


26.9 


0.21 


41.6 


2,4 


1306-09 


0.464 


27.1 


2.72 




1 


1328+307 (3C 286) 


0.849 


28.2 


22.0 


44.3 


2,3 


1345+125 (4C 12.50) 


0.122 


26.0 


0.16 


43.0 


2,3 


1358+624 (4C 62.22) 


0.429 


26.9 


0.35 


41.8 


2 


1443+77 (3C 303.1) 


0.267 


26.0 


6.22 


42.9 


2,3 


1634+628 (3C 343) 


0.988 


27.7 


1.45 




2 


1637+626 (3C 343.1) 


0.750 


27.4 


1.57 


43.3 


2,3 


1807+698 (3C 371) 


0.050 


25.0 


~ 1 


40.7 


3,6 


1934-63 


0.183 


26.7 


0.17 




1 


2342+821 


0.735 


27.4 


1.17 


43.1 


2,4 


2352+495 


0.237 


26.3 


0.17 


41.3 


2,4 



Note. -All quantities are caluculated assuming Hq — 75 kms 1 Mpc 1 and qo — 0. Col. (1): Object name. Col. (2): rcdshift. Col. (3): log of radio 
power at 5 GHz. Col. (4) linear size of the radio jets. Col. (5): log of [OIII]A5007 luminosity. 

Referenccs.-(l) Tadhunter ct al. 1993, (2) O'Dca 1998, (3) Gcldcrman & Whittle 1994, (4) Lawrence et al. 1996, (5) Kawakatu, Nagai & Kino 2008, 
(6) Lister ct al. 2001, 
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Table 2. Emission-line ratios relative to H/3 flux 



Name 


[Oil] 


[NelllJ 


[Olll] 


[Olllj 


lull 


[hill] 


[Sllj 


References 




A3727 


A3869 


A4363 


A5007 


A6300 


A6716 


A6731 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


0023-26 


4.90 


0.51 


0.10 


3.40 


2.70 


2.48 


3.10 


1 


0316+413 


2.95 


0.92 


0.38 


10.9 


6.61 


3.37 


2.67 


3 


0345+337 








8.33 


2.00 






2 


0404+769 




0.33 




3.66 








3 


0518+165 








0.93 








2 


0605+480 


5.25 






4.74 


2.26 


3.17 


2.71 


3 


0710+439 


2.54 


0.44 


0.28 


6.10 


1.33 






3 


0954+658 


0.45 


0.47 


0.33 


7.48 








3 


1031+567 


3.16 


0.98 




7.20 


5.80 






3 


1306-09 


4.64 


0.60 


0.06 


3.57 








1 


1328+307 


0.48 


0.57 


0.39 


4.34 








2 


1345+125 


2.19 


1.06 


0.38 


6.25 


0.94 






2 


1358+624 


2.49 


0.78 


0.10 


8.69 


5.61 


3.49 


5.04 


3 


1443+77 


4.58 


0.50 




9.00 








4 


1634+628 


2.06 


1.08 


0.18 










3 


1637+626 








6.33 








2 


1807+698 


1.10 






3.66 


2.34 


1.06 


2.01 


3 


1934-63 


3.12 


1.44 


0.26 


7.21 








1 


2342+821 


0.94 


0.59 




9.98 








3 


2352+495 


3.21 


0.80 




8.56 


3.88 


3.25 


3.01 


3 



Note. -Col. (1): Object name. Cols. (2)-(8): Emission line ratios relative to H/3 flux. 

Referenccs.-(l) Morganiti et al. 1997, (2) Gcldcrman & Whittle 1994, (3) Lawrence et al. 1996, (4) Labiano ct al. 2005 
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logio([OII]X3727/[OIII]X5007) 

Fig. 1. — Diagram of [Ol] A6300/[Om] A5007 vs. [Oll]A3727/[Om]A5007. The observed line ratios from the literature are presented for 
powerful YRL AGNs (red circles) and radio-quiet Seyfert 2 galaxies (black crosses). Solid and dashed lines represent our photoionization 
model calculations using the SED with BBB and the SED without BBB, respectively. 
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Fig. 2. — Frequency distributions of emission line flux ratios [Ol]A6300/[Oni] A5007, for powerful YRL AGNs (red) and Seyfert 2 galaxies 
(blue). 
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Fig. 4. — The optical emission line ratios ([Oi] A6300/[Om] A5007 and [Oii]A3727/[Oin]A5007) as a function of the power-law index of 
RIAF-like SEDs, opl (left), and the characteristic temperature, Tbb (right) for a standard accretion disk, assuming rtjj = 10 4 cm - 3 . The 
solid lines represent the case of log U = —3.0, while the dotted lines denote that of log U = —3.5. The blue (red) solid line represents the mean 
value of log([Ol]A6300/[Oni] A5007) for powerful YRL AGNs (Seyfert 2 galaxies). The blue (red) dotted lines correspond to the dispersion of 
log ( [Oi] A6300/ [Oni] A5007) . 
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Fig. 5.— Diagram of [Om] A5007/[Oin] A4363 vs. [On]A3727/[Om]A5007. Symbols and lines are the same as in Figure 1. 
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Fig. 6. — Frequency distributions of the emission line flux ratios of [Om] A5007/[Oni] A4363, for powerful YRL AGNs (red) and Seyfert 2 
galaxies (blue). 
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Fig. 7.— (a) Diagram of [Oi] A6300/[Om] A5007 vs. [On] A3727/[Om]A5007, (b) diagram of [Oi] A6300/[Om]A5007 vs. 
[Sii]A6717+6731/[Oiii]A5007, and (c) diagram of [Ol]A6300/[Oni] A5007 vs. [Nelll]A3869/[Olll] A3727. Predictions of the shock model (Dopita 
& Sutherland 1995) are represented by blue lines. The blue solid lines denote the pure-shock model (150km/s ^shock 

< 500km/s), while 

the blue dotted lines denote the models considering the effect of precursor (200km/s < i> s hock < 500km/s). Chemical composition of solar 
abundances and the magnetic parameter of B/^/n = 2/iGcm 2 / 3 are assumed in these models (see Table 1 in Dopita & Sutherland 1995 for 
details). Symbols and line are the same as Figure 1. 
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Fig. 8. — Schematic sketches for the properties of accretion disk and the obscuring materials of the RQ Seyfert 2 galaxies (upper) and the 
powerful YRL AGNs, i.e., CSOs and MSOs (lower). The light gray circles denote the narrow line clouds emitting [Oin]A5007 and the dark 
gray circles represent the narrow line clouds emitting [Oin]A4363. 



